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Abstract 

In this paper we study the nonabehan, gauge invariant and asymptotically free quantum 
gauge theory with a mass parameter introduced in [8]. We develop the Feynman diagram 
technique, calculate the mass and coupling constant renormalizations and the effective 
action at the one-loop order. Using the BRST technique we also prove that the theory is 
renormalizable within the dimensional regularization framework. 
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1 Introduction 

In paper [8] we introduced a new classical gauge invariant nonlocal lagrangian generating a local 
quantum field theory which is reduced to several copies of the massive vector field when the 
coupling constant vanishes. This unexpected phenomenon is similar to that for the Faddeev- 
Popov determinant in case of the quantized Yang-Mills field. Recall that being a priori nonlocal 
quantity the Faddeev-Popov determinant can be made local by introducing additional anticom- 
muting ghost fields and applying a formula for Gaussian integrals over Grassmann variables. 
Similarly, in case of the lagrangian suggested in [8] one can introduce extra ghost fields and 
make the expression for the generating function of the Green functions local using tricks with 
Gaussian integrals. 

In this paper we develop the Feynman diagram technique, calculate the mass and coupling con- 
stant renormalizations and the effective action at the one-loop order for the theory introduced 
in [8]. Using the BRST technique we also prove that the theory is renormalizable within the 
dimensional regularization framework. 



2 Preliminaries 

First we fix the notation used throughout this paper. Let G be a compact simple Lie group, q 
its Lie algebra with the commutator denoted by [■,•]. Denote by tr the Killing form on g with 
the opposite sign. Recall that tr is a nondegenerate invariant under the adjoint action scalar 
product on g. Let a = 1, . . . ,dimg be a linear basis of q normalized in such a way that 
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tj:{t°',t^) = 5°"^. Note that for technical reasons the normalizations of scalar products in this 
paper are slightly different from [8] and from the standard ones (see [5]). 
Since the tensor formulas that will appear in this paper are even more complicated than in 
case of the Yang-Mills field we shall frequently use the invariant tensor notation and a shorten 
notation for scalar products of tensors. Denote by 0*^(0) the space of fl-valued differential 
k- forms on the standard Minkowski space equipped with the metric g^^, 500 = 1, ffii = — 1 for 
i = 1,2, 3, and gij = for i 7^ j. 

We shall use the scalar product < •, • > on the space fl''{g) defined by 



< U!,u)' >= J tr(wA,*a;'), lo,u)' G ^''{q). 



(1) 



where * stands for the Hodge star operation associated to the metric on the Minkowski space, 
and we evaluate the Killing form on the values of wi and *U2 and also take their exterior 
product. 

Let A e i^^io) be the fl-valued gauge field (connection one-form on the Minkowski space), 



A = Af.dx'' = Al^dx", 
and F G ^^{q) the strength tensor (curvature) of A, 

F = dA-^[AA,A]. 

Here g is a coupling constant and as usual we denote by [^A,j4] the operation which takes 
the exterior product of g-valued 1-forms and the commutator of their values in fl. In terms of 
components we have 

F = F^^dx" A dx"" = ^(9^A, - d^A^ - g[A^, A,])dx^' A dx\ 

We shall also need the covariant derivative operator dA ■ i^*^(0) ^''^^{q) associated to the 
connection A, 

dAijJ = duj — g[A/\,ui\. 

The operator d*^ adjoint to dA with respect to scalar product (1) is equal to (— l)'^(''~'^)+i 

The gauge group of G- valued functions g{x) defined on the Minkowski space acts on the gauge 

field A by 

A^-dgg-^+gAg-\ (2) 
g 

This action generates transformation laws for the covariant derivative and the strength tensor, 

dA^gdAg~^, (3) 
F ^ gFg-\ (4) 

In formula (3) we assume that the gauge group acts on tensor fields according to the repre- 
sentation of the group G induced by that of the Lie algebra q, and the r.h.s. of (3) should be 

regarded as the composition of operators. 

We shall also need a covariant rough d'Alambcrt operator Ua associated to the gauge field A, 

Ua = D^D'', D^ = d^-gA^, 

Here and thereafter wc use the standard convention about summations and lowering tensor 
indexes with the help of the standard metric on the Minkowski space. 
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The covariant d'Alambert operator can be applied to any tensor field defined on the Minkowski 

space and taking values in a representation space of the Lie algebra g, the g-valued gauge field 
A acts on the tensor field according to that representation. Note that the operator Da is scalar, 
i.e. it does not change types of tensors. 

Formula (3) implies that the covariant d'Alambert operator is transformed under gauge action 
(2) as follows 

□a ^ gOAg-^. (5) 

For the zero connection ^ = we simply write do = d and Do = □, the usual d'Alambert 

operator. 

Using the tensor notation the "localized" action for the massive gauge field A introduced in [8] 
can defined by the formula 

S' = ~<F,F>+j^<aA^,^>~ <^,F>+ (6) 

3 

+ ^ < dAV^dAVi > ■ 

i=l 

Here $ € ^^{q) is a bosonic ghost field in the adjoint representation of g; r]i,rji, i = 1, 2. 3 are 
pairs of anticommuting scalar ghost fields in the adjoint representation of g. In formula (6) g 
should be regarded as a coupling constant and m is a mass parameter. From (3), (4) and (5) 
it follows that action (6) is invariant under gauge transformations (2). 

In formula (6) we use normalizations slightly different from those introduced in [8] . For instance, 
the coefficient in front of the last term in (6) is different from that which we used in [8] . But 
this coefficient is not important for the study of renormalizability of the theory. 
To define the Green functions corresponding to the gauge invariant action S' we have to add 
to action (6) a gauge fixing term Sgf. As it was observed in [8] the most convinient choice of 
the gauge fixing term is 

A TYX^ 

Sgf = -- < d*A,d*A > — — < D~'^d*A,d*A > + < djr],dAr] >, 

where r],!] is a pair of anticommuting scalar ghost fields in the adjoint representation of g, and 
we introduced dimensionless parameter A to consider different gauge fixing conditions (in [8] 
we only discussed the case when A = 1). In this paper we shall also assume that all inverse 
d'Alambert operators are of Feynman type. 

The total action S = S' + Sgf that should be used in the definition of the Green functions takes 
the form 

1 1 Tfl 

s = -- < F,F > +— <nA<^,^ > - J < '^,F > + < dr],dAV > + (7) 

3 A 
+ ^ < dAfii,dAVi > ~2 ^ d*A,d*A > — — < n~^d*A,d*A > . 

i=l 

Now consider the expression for the generating function Z{J) of the Green functions for the 
theory with action (7) via a Feynman path integral, 

3 

Z{J) = / V{A)V{<^)V{n)V{l]) [] V{r^,)V{r^,) exp{z(5+ <J,A >)}, (8) 

■' i=l 

where J is the source for A. Here and thereafter we assume that the measure in Feynman 
path integrals is suitably normalized. We also assume that the Feynman path integral over $ 
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in (8) is only taken with respect to the hnearly independent components ^^v, H < f of the 
skcw-symnictric tensor 

Observe that in the r.h.s. of formula (8) all the integrals over the ghost fields are Gaussian. 
The Gaussian integrals can be explicitly evaluated (see [8] for details). This yields 

/I ^ 
P(A)det(d*d^|oo(g))exp{i(-- <F,F>-'^< n-/F,F> + (9) 

+ <J,A> --< d*A,<fA > < n-^d*A,d*A >)}. 

The r.h.s. of (9) looks like the generating function of the Green functions for the Yang-Mills 

theory with an extra nonlocal term, < O'^^F^F >, and in a generalized gauge (see [3], 

Ch 3, §3). Although the action Sm that appears in formula (9), 

Sm = -l<F,F>-^ < D^'F,F> ~ < d*A,d*A > - (10) 

2 

< a-^d*A,d*A >, 

2 

is not local the generating function Z{J) for this action is equal to that for action (7). It is the 
phenomenon that was observed in [8]. 

The purpose of this paper is to study generating function (8) in detail. In the next section we 
shall discuss the Feynman diagram technique for action (7). 



3 Feynman diagram technique 

In order to develop perturbation theory for generating function (8) we have to introduce addi- 
tional sources / € ^^{q), d « = 1,2,3 for the ghost fields 77, rj, rji, rj^, i = 1,2,3, 
respectively. We also need the full generating function Z{J, I, ^, ^, f*) for the Green functions, 

_ _ f 3 
Z{J,I,^,C,e,^')= / V{A)V{<^)V{n)V{rj)Y[V{rh)V{r],)x 

i=l 

X exp{i(S'+ < J, ^ > +i </,$> + < > + <■?,??> + (11) 

3 

+ Y.^<t,rii> + <e,^>))}- 

i=l 

Recall that the connected Green functions are defined as the coefficients of the generating series 
G{J,I,^,l,€,t) related to Z(J,/,^,|,e,f ) by 

Z{.J, I, e, I, e,t) = exp(G( J, I, ^, I, f , ?)). 

We shall first calculate the free propagators for the theory with generating function (8). The free 
propagators arc the connected Green functions corresponding to the free generating function 

- ■ r ^ 

Zo{J,I,^,te,f) = I V{A)V{^)V{rj)V{r])l[V{r]i)Vif],) x 

x exp{i(S'o-|- < J,A> I,^ > + <^,T]> + <^,r}> + (12) 

3 

+ J2{<t,Vi> + <e,rii>))}, 

i=l 
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where 



1 1 fn 

So = --< dA,dA > +— < > -— < ^,dA> + < djj,dr] > + (13) 

^ A 
+ ^ < df]i,dr]i > -- < d*A,d*A > — — < n-'^d*A,d*A > 

is the unperturbed action corresponding to (7). 

Feynman path integral (12) is Gaussian, and we can exphcitly evaluate it. First we integrate 
over the ghost fields in (12). This yields 



Zo{J,I,^,^,e,f) = j V{A)e^v{i{-\ < {d*d+dd*)A,A> 



- ^^ < dd*A,A> < {d*d + dd*)D-'^A,A> - < mn-'^dA,I> + (14) 
1 - ^ • 

+ <j,A>--< D-'i,i> + < n-'^,^>+J2 < n-T,e >)}• 

i=l 

Now using the Hodge formula, dd* + d*d = —D, in the exponent in (14) and integrating over 
A we finally obtain 

Zo{J, I, e,?, C,t) = exp{i(-i < (□ + + (1 - X)dd*)-^J, J>- 

< (D-i + m^n-^dd*{n + + (1 - X)dd*)-'^)1, 7 > - (15) 
- < m(n + + (1 - X)dd*)-'^n-^d*I, J> + 

3 

+ <n-i?,^>+^<n-T,f >)}, 

where all the inverse d' Alambert operators are of Feynman type. We see that the free propagator 

for the gauge field A coincides with that for a massive field. When A = 1 it is reduced to 
the propagator for the free massive vector field, i{D + m'^)~^ . The fact crucial for the mass 
generation is that free action (13) contains a cross term, — ^ < ^,dA>, which generates mass 
in the free propagator of the gauge field. Note that the free propagator is also not diagonal. It 
contains the cross term < m,{D + + (1 — X)dd*)~^D~^d* I, J >. 

We shall develop perturbation theory in the momentum space for the Fourier transforms of the 
Green functions. We normalize the Fourier transform as follows 



{2n) 

The Fourier transforms of the propagators arising from (15), along with the corresponding 
diagrammatic symbols, are summarized in the following table: 




i(n + m2 + (l-A)rfd*)-^ = 



ab 



9^v 



+ ■ 



(l-A-i)fc^fc, 



— m? + ie [k 
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A* < 7 k V <5 
It^f > i{n-^ + m'^n-^dd*{n + + (l - X)dd*)-'^) = 



a 

-US"'' r 
fc^ + is 



9ij.i'9i5 + _|_ j£)(yt2 _ ^2 _|_ jg.-) idni'kjks Qjuki^kg Qnskjk^ + Qjskfik^) 



k ^ <1 



/^^ im(n + m2 + (l-A)drf*)-in-irf* = 



(16) 



2m5 



ah 



a 



1« 



A;2 + i. 



In the diagrams above we indicate the sources coupled to the corresponding propagators in 
formula (15). We assume that in coordinate formulas containing summations over indexes of 
skew-symmetric (2,0)-type tensors all summations are only taken over the indexes correspond- 
ing to the independent components of the skew-symmetric tensors; for instance, in formulas 
containing summations over the indexes of the source Z^" we take sums over jj, <u. This gives 
some extra numerical factors in the formulas for the propagators. For instance, we have 

<ui,J>=A I ^ir{uj^^,uj'i^'')d^x, £n^{Q). 
■' iJ.<f 

In the formulas for the Fourier transforms of the propagators given above we also omit, as usual, 
the delta functions of the total external momenta. 

The perturbation theory for generating series (8) is based on the following formula 



ZiJ) = 



(17) 



I vf — — — — — — 

^""P I [i6j' isr ise i6t iS^' ' isf 



/=^=|=|i=^=0'^0('^)-f)C)^)C)^ )' 



where 



V{A, r], r], r]i, r?j) = S{A, rj, r], r]i,r]i) - So{A, r?, r], r]i,r]i) 



is the perturbation of the free action 5*0, and in formula (17) all variational derivatives with 
respect to the Grassmann variables ^, S^,- (^, ^i) arc right (left), respectively. 
To complete our study of the perturbation theory based on formula (17) we have to calculate the 
contributions of the vertices corresponding to the terms of the perturbation iV. These terms and 
their contributions evaluated in the momentum representation, along with the corresponding 
diagrammatic symbols, are listed below. 
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< dA, [AA,A] >, 
gC'""={2n)^Sip + q + r){g^.{p - q)p + g,p{q - r)^+ 




-^<[AA,A],[AA,^]>, 
^ -ig2(27r)45(p +q + r + s){C'^''C''''{g^pg^^ - g^ag.p)+ 



p V c 



a < /? 




%<d^^,[A^,^>, 



|(27r)^C"'"=5(p + ? + r)(p^ - q^,){gai9p5 - gMQas) 



7 < (5 



q 



a 



g<[^^,$],[A^$]>, 



r 



s d 



x(5'a75/3<5 - 9a590-i) 
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^<$,[AA,A]>, 



(18) 



-ig<dfj, [A,T]] >, 

'>%2Tr)^6{k+p-q)p^ 



-2ig<dfji,[A,rii] >, 
-2gC'''"={2T,Y5{k+p-q)p^ 



ig^ < [A,T)i]h,[A,r]i] >, 



2ig2(27r)*5(p + q + r + s)g^^{C''"-''C''^'^ + C'bcQead 



(19) 



" d 

In the coordinate formulas above C'^^'^ are the structure constants of the Lie algebra g, t^] = 
(jabc^c^ In the diagrams containing ghost lines we assume that the variables rj, rji correspond to 
the outgoing ghost lines. All the coordinate formulas include the corresponding combinatorial 
factors. Each diagram constructed with the help of vertices and lines given above gives a 
contribution to the generating function Z{J). To calculate the contribution of a connected 
diagram one has to integrate the expression corresponding to the diagram with the measure 
over all internal momenta and multiply the result by ^ g' , where I in the number of the 
T]- and 77j-ghost loops and S is the order of the symmetry group of the diagram. 
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4 Effective action at the one-loop order and one-loop renor- 
malization 

In order to calculate the Green functions at the one-loop order we shall use a formula expressing 
the effective action at the one-loop order via a Gaussian Feynman path integral (see [5], Sect. 
12.3). We recall that the effective action T{A) is defined as the Legendre transform of the 
generating function G(J) = In Z{J), 

iT{A) = G{J)-i<J,A>, A=^. 

loJ 

At the lowest order the effective action coincides with the classical action of the theory. The 
first quantum one- loop correction T^{A) can be represented as follows 

iV^{A) = (20) 
= \n\^jv{A')V{^')V{r)')V{rj')f{v^^^^^ 

where S^{A] A' ,<^' ,r}' ,r]' ,ri^,%) is the quadratic term in the Taylor expansion of the classical 
action S{A + A', r]',r]', rjljTj^), with respect to the primed variables at point {A, 0, 0, 0, 0, 0). 
Prom formula (7) using the definition of the curvature we deduce that 

52 = I < *[*F(^)A, A'], A' >-^<dA* dAA', A'> + 

1 7TI 

+ — < D^*', $'>--<$', dAA' > + <rf, d*dAv' > - (21) 

-y2<v'i,OAv'i >-^< dd*A',A' > < dn-^d*A',A' > . 
i=i ^ ^ 

Now Gaussian integrals in formula (20) can be easily evaluated. This leads to the following 

expression for r^(A) 

ir\A) = ln{det(-n-^d*dA)|ao(0) (det [{a + m'^ + {1 - X)dd*)-'^x (22) 

x(g * [*F{A)A, ■ ] - d*AdA - m''d*An-/dA - Xdd* - m^dn-irf*)|ai(0)])"h. 

As usual wc normalize the one- loop effective action in such a way that T^{A) = when g = 0. 
Recalling that In det X = trlnX for any operator X we can rewrite formula (22) in a form 
suitable for calculations, 

iT^{A) = tr{ln(-n-id*rfA)|ao(g) - i In [(□ + m'^ + {l- X)dd*)-^ x (23) 
x(g * [*F{A)A, ■ ] - d*AdA - m^d^n^^^dA - Xdd* - m'^dD-^d*)\n^g)\ }. 

Formula (23) is the most explicit general expression for the effective action at the one-loop 
order. The first term in the r.h.s. of (23) is the contribution related to the Faddeev-Popov 
determinant. 

Now using formula (23) we shall calculate the one-loop correction r2 to the two-point function 
for the gauge field A. In practical calculations we shall use the momentum representation. 
First we apply the Wick rotation to the time components of all momenta and of the gauge field 
entering formula (23), 

k°^ik°, A°{k)^iA°{k). 
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We shall use dimensional regularization to preserve gauge invariance, so that all the integrals 
with respect to momenta are taken over M'' with respect to the measure • In the regularized 

expressions evaluated in K'^ the coupling constant g should be replaced with gfi^ ,e = A — d, 

where jjL is an arbitrary mass scale, so that g remains dimensionless in R'^. 

The contribution {r\)ghost to the two-point function related to the Faddeev-Popov determinant 

tr{ln(-n-id*dA)|no(g)} = tr{ln(7 - gn-^d^[A^, ■ ])|oo(g)} 

in the r.h.s. of (23) is standard and coincides with that for the pure Yang-Mills field. Indeed, 
expanding the r.h.s. of the last formula with respect to g we obtain 

i{T\)ghost = -y tr {n-^d^ o [A^ . ] o D-'d, o [A'', ■ ]|oo(g)) , (24) 

where the symbol o stands for the composition of operators. The trace in formula (24) can be 
evaluated using the momentum representation, with the conventions about the Wick rotation 
and the dimensional regularization introduced above, 

Hi2W- 2 (p + 5)V (27r)<^(27r)<i' 

where • denotes the standard scalar product of vectors in d-dimensional Euclidean space W^, 
= q ■ q, etc., and the constant C is defined with the help of the structure constants of the 
Lie algebra g, 

^ja6 ^ ^ ^acd^bcd 

c.d 

The contribution {TI)a to the two-point function related to the second term 

tr{-i ln[(n + m'^ + {l- A)rfrf*)-^(g * [*F{A)A, ■ ] - d^d^ - (26) 

-m'^d*An;^^dA - Xdd* - m'^dn-^d*)\n^^)]} 

in the r.h.s. of (23) can be calculated in a similar way. For simplicity we shall only consider 
the case when A = 1. Expanding the r.h.s. of (26) with respect to g and using the resolvent 
formula for the d'Alambert operator D^^, 

= (/ - gD-\d^ o [A^ ■ ] + [A^ • ] o 5.) + g^n-i o [A^ . ] o [A^, ■ ])-'a-\ 

we deduce that 

z{Tl)A=^tv{^X',-X2)\nH,), (27) 
where the operators Xi and X2 are defined by 

Xi = (a + m^)-! {d* o [^A, • ] + [AA, ■]*od- 
-m^d*D-'\d^o[A'',-] + [A'',-]od^)D-^d + m^[AA,-]* oD-^d+ (28) 

+in^d*a-^ o [AA, ■ ] + *[*idA)A, ■ ]) , 

X2 = {D + m'^)-\-[AA, ■ ]* o [AA, ■ ] + m^d*^-'^ o [A'^, • ] o [^^, • ] o Q-irf - 

-m''d*n-\d^ o [A'', ■ ] + [A^ • ] o d^)D-\d, o [A-', ■ ] + [A", ■ ] o d,)a-^d + 

+m^[AA, ■ ]* o □-1(9^ o [A^', ■ ] + [A", ■ ] o d^)a~^d + (29) 
+m^d*a-\d^ o [A^, ■ ] + [A-", ■ ] o 9,)n-i o [AA, ■ ] - 

-m^[AA, • ]* o o [AA, • ] - i * [*[AA, A], • ]). 
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In formulas (28), (29) the superscript * for operators always denotes the operator adjoint with 
respect to the scalar product on il^(g). 

Now we have to calculate the traces of the operators and X2 using the dimensional regu- 
larization as above. After tedious algebra we get 



4 {p + 2g) • (p + 2q) ■ A-{-p){{d - 2){{p + q) ■ qf + {p + qfq^) 



(P + 2g) ■ A-{p){id - 2)q ■ A-{-p) {p + q)-q- q^p + q) ■ A-{-p)) ^ 
g4(p + g)2((p + q)2 +^2) 

+2 ^4((p+^)r+^2) ((^ - 2)9 • q ■ + q'A'^i-p) ■ - (30) 

2 (p + 2q) ■ l°(p)(g • >(-p)p^ -p-qp- 1°(-P)) , 



+4 



ip + qYq^{q^ + 'rn?){(p + qY + w^) 



q^{{p + qf +w?) 



p ■ A'^ip) p ■ - p'A'^ip) ■ 

trx,=c(d-i) / t ' ^ • - (31) 

^ ^' (27r)'^ (27r)'' ^ {p + qfq^{q^ + tv?) 



^^ p-A-j-p) i2q+p)-A-ip) ^ 
{p + qYq^{q^ + m^) 



As a function of e, has a simple pole singularity when e = 4 — d ^ 0. To find the counterterm 
corresponding to the function r2 we have to calculate its residue at point e = using formulas 

(25), (27), (30) and (31). 

Note that we use the minimal subtraction procedure to regularize the efi'ective action. Therefore 
the divergent part of r2 , that we use to find the counterterms, is the principal part of the Laurent 

series of r2(e) at point e = 0. 

The integrals over q in (25), (30) and (31) can be evaluated with the help of the following 
formulas 



1 



(27r)'* {{p - qf + mf )"(g2 + mlY 

d\ ,-1 



I 

..\^^rArr f daa^-\l - af-\a{l a)p' + ami + (1 - a)ml)i-^-^, 

(477)2 r(n)r(p) Jq 

[ dS qj. ^ 

J (27r)d((p-g)2+m2)"(g2 + rn2)P 

(47r)2 r(n)r(p) Jq 



(32) 



(33) 
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r cfiq 
J J2^W^' 



(34) 



(47r)i' r(n)r(p) Jo 

+ (l-a)TO^)^-"->^p^ + 

r(n)r(p) Jo 

+(l-aK)^+i— ^-^i, 

where F is the Euler gamma function. 

Finally recalling that for e ^ the principal part of the Laurent series of the function F(— n+e), 
n = 0,1, 2,... is equal to we obtain an expression for the divergent part {T2)div of the 

function 

Returning to the Minkowski space by the inverse Wick rotation we can rewrite the last formula 
in the configuration space as follows 

(^2)div = T^^il < dA,dA > +2m2 < A,A>). (35) 
(47r)^e 3 

The divergent term (35) can be eliminated by adding a counterterm ^2 to nonlocal action (10), 

Sl = --^^{1 <dA,dA> +2m^ <A,A>). (36) 
(47r)^£ 3 

This is, of course, equivalent to a modification of the coefficients of local action (7). Actually 
only the coefficients in the expression for the free action (13) have to be renormalized, and the 
expression for the renormalized free action (6*0)2 takes the form 



(50)2 = ~Za < dA, dA > +^ <□$,$> -^Zm <^,dA> + <d^,dr]> + (37) 

^ 1 
+ Y]< d%,drn > -7: < d*A,d*A > -^Z^ < n^U*A,d*A >, 

where 



1 

By rescaling field A in formula (37), A° = Z^A, we obtain 

{So)l = -\< dA°,dA° > +^ <□$,$> -jZmZ^^ <^,dA° > + <an,dr]> + (39) 
+ ^ < dni,dr]i > < d*A°,d*A° > -I^Z^Z^i < a-^d*A°,d*A° > . 
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Now Prom, formulas (38) and (39) we can find the mass renormalization, 



mo = mZ^Zj^ ^ = 



m 1 - 



11^2 



) 



(40) 



167r2 6 e 



where mo is the bare mass. 

Note that since the possible renormalization constant of the ghost field $ is canceled in the 
Gaussian integral in formula (12) the renormalization of the wave function of $ is not required 
for calculating the mass renormalization. 

To calculate the one-loop coupling constant renormalization one has to find the divergent part 
of the three-point function T\ at the one-loop order. Looking at formulas (28) and (29) one 
should expect that technically this problem is extremely complicated. 

In the next section we shall prove the universality of the coupling constant renormalization and 
Slavnov- Taylor identities for the theory with generating function (8). Using the universality of 
the coupling constant renormalization one can also find the one-loop coupling constant renor- 
malization with the help of the one-loop corrections to the three-point function represented by 
diagram (19). This calculation is standard and is completely similar to the case of the pure 
Yang-Mills field (see [3], Ch. 4, §1). We only note that in case of the dimensional regularization 
the mass in the propagator for the gauge field does not add new divergent terms to the quantum 
corrections. The coupling constant renormalization computed in this way coincides with that 
for the pure Yang-Mills field, 



where go is the bare coupling constant. 

One can also observe that the coupling constant renormalization in the theory with generating 
function (8) calculated within the dimentional regularization framework coincides with the 
coupling constant renormalization for the pure Yang-Mills theory to all orders of perturbation 
theory. Indeed, by the results of [4] in the case of dimensional regularization the coupling 
constant renormalization is independent of the mass for dimensional reasons. Therefore the 
coupling constant renormalization in the theory with generating function (8) is the same as in 
the massless case, i.e. in case of the Yang-Mills field. 

From formulas (38), (40) and (41) one can calculate the corresponding renormalization group 
coefficients at the one-loop order. An elementary renormalization group analysis of the theory 
with generating function (8) can be found in [8]. 

5 Renormalization 

To complete our investigation of the theory with generating function (8) we have to prove its 

rcnormalizability to all orders of perturbation theory. 

By simple dimensional counting all the integrals that enter the perturbative expressions for the 
Green functions can be dimensionally regularized. Indeed, the degree of divergence uj{G) of 
any diagram G that appears in the perturbative expansion for generating function (8) can be 
computed as follows: 



where E is the number of the external lines in the diagram G, 5 is the total power of external 
momenta factorized from the corresponding integral, and the sum in the last term is taken over 
parts (16) of the free propagator which appear as internal lines in G and over vertexcs (18). 
From formula (42) we deduce that divergences may only appear in the Green functions with 
E = 2,i, 4, and the theory is renormalizable. 




(41) 




(42) 
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We have to prove yet that the gauge invariance of the theory is preserved by the renormalization. 
As usual this is achieved with the help of the corresponding Slavnov- Taylor identities that we 
are going to derive now. 

The proof of renormalizability of the theory with generating function (8) is completely analogous 
to that for the Yang-Mills field (see [5], Sect. 12.4). In the proof we shall use the BRST 
technique. We shall consider the case of arbitrary g-valued gauge fixing condition J^{A) linear 
in A, J^{A) = (j^ab-^ix- ^^^^ situation is slightly more general than the one discussed above for 
action (7), and the new total action including the gauge fixing term takes the form 

1 1 Tfl 

Sy. = --<F,F>+—< > -_ < $,F > - < 7],Mv > + (43) 

3 A 
+ ^ < dAr]i,dAVi > - 2 < >> 

i=l 

where M is the Faddeev-Popov operator for the gauge fixing condition J^{A). 

The BRST transformation s can be defined similarly to the case of the Yang-Mills field ([5], 

Sect. 12.4.2), 

sA = dATj, srj = XJ^{A), 

s$ = g[r], $] , sr]i = g[r], rji] , (44) 

srj = ~^C°-'"'r]''T]'', sr]^ = g[v, Vil- 
li we introduce the ghost number gh for the fields as follows 

ghA = ghi> = gh?7i = ghrji = 0, 

gh?7 = -1, ghr? = 1, 

then s becomes a right superderivation of the algebra of the fields, and 

s'^A = = s^?7 = s'^rii = 3^% = 0, s^l) = 0. 
We also define three ghost numbers gh^, i = 1,2,3 by 

gh'iA = gh.$ = gh-r? = gh-^ = 0, 

We shall find the Slavnov- Taylor identities for generating function G(J, /, S,\k,L,N,P^,7^), 

_ . _. r 3 

exp G{J, /, e, ^, e, f , K, L,N,P\P')= / V{A)V{^)V{r,)V{rj) JJ -DiViMv^) x 

i=l 

X exp{i{Sj^+ < J,A> I,^ > + <'^,r] > + <^,r] > + (45) 

3 _. I 

+ ^(< ^\vi > + < CVi >)+ < K,sA> - < L,ST]> +-< N,s'i> > + 

i=l 

3 

+ ^(< P\sT]i > + <P\sr)i >)}, 
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where K,L,N,P^,P arc the sources coupled to sA, srj, s$, srji and srj^, respectively. The 

generating function G( J, /, ^, ^, C , K, L, N, P% P*) is reduced to G( J) when all the sources, 
except for J, are equal to zero. The new sources have the following ghost numbers 

ghK = ghN = ghP* = ghP' = 1, ghi = 2, 

gh^i^ = gW^N = gW.L = 0, gh^P^ = Sij, gh^F = -6ij. 
We shall also need the dimensions of the fields and of the sources, 

dimj4 = dim$ = dimr/ = dimry = dimr^i = dimry^ = 1, 

dimJsT = dimL = dimiV = dimP, = dimPj = 2, 

dim J = dim/ = dim^ = dim^ = dim^' = dim^ = 3, 

the dimension of the mass m being equal to 1, and the coupling constant g is dimensionless. 

To derive the Slavnov- Taylor identities we have to apply the BRST transformation (44) to the 
variables of integration in the r.h.s. of formula (45). Observe that both action (43) and the 
measure in the Feynman path integral (43) are BRST-invariant. The former statement can 
be checked directly, and the latter one is true since the BRST transformation is nilpotent. 
Therefore the BRST coordinate change in (43) yields the following identity 



/ V{A)V{^)V{n)V{r])f[V{r]i)Virji){< J,sA> +^ < I,s^ > + < ^sv > - 

3 _. 1 

- < sT],^ > +^{<t ,sr]i > + < srii, C >)) exp{i{Sj^+ < J,A> +- < I,^ > + 

i=l 

3 

+ <l,V> + <^,r]>+ X](< Vi> + <C,Vi>)+ < K,sA> - < L,sr]> + 

i=l 

1 ^ -• 

+-< N,s^> + Y,{< P\ sr]i> + < P\ srJi >)} = 0. 



i=l 

This identity can be rewritten as 

(<^'^>+4<^'^>-<^'I^>-^<^'^(z^)> + 

+ E(< ^ > - < ^ >)) exp G{J, I, C,t, K, L, N, P^ p') = 0, 

i=i i6P 

which is, in turn, equivalent to 

+ E(< ^ > - < e, ^ >))G{j, I, a, e,T, K, L, N, p\T) = 

due to linearity of T. 
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The Slavnov- Taylor identities (46) are supplemented with another identity following from the 
fact that the Feynman path integral in the r.h.s. of formula (45) is invariant under translations 
of the ghost variable of integration, 

r] i-^r] + Sfj. 



Prom this fact we infer that 

3 



/ V{A)V{<^)V{r,)V{r}) \{ V{rh)V{r}i){-Mii + exp{i(5^+ < J, ^ > +1 < /, $ > + 

i=i 

3 

+ kI,^) > + < i,r) > +X](< m> + < C,Vi >)+ <K,sA> - <L,sr]> + 

i=l 

1 -• 

+ - < iV, s$ > + ^(< P^ sr]i> + < P\ srJi >)} = 0, 



4 



or, since M.r] = sJ^ = (t^abi^^'n ■~ ^ ' 



(C - • T^) exp G{J, I, ^, K, L, N, P\ p') = 0. 

For connected functions the last identity gives 

• ^' ^' e,t,K,L,N,P\T)= e (47) 

In the proof of renormalizability it is more convinient to use the proper functions instead of 
the Green functions, and instead of the generating function G{J,I,^,^,^'^,^,K,L,N,P^,^) 
we shall consider its Legendre transform T{A, rj, rj, rji,rj^, K, L, N, P', P*), 

r(^, ^,7i,rj,v,,rj„ K, L, N, P\p') = -iG{J, I, ^, C, K, L, N, P\p') - (48) 

1 - ^ -■ 

-<J,A>+-<I,^> + <^,r] > + +X](< > + < ^''^i >)' 

where 

, 5G SG _ 6G ^ SG SG _ SG 

10 J iS^ «o-/ 

and the variables if, L, N, P', P' are passive. 
Since 

^ Srj 

and 



(5r 

^~ 


I = - 


ST -t 


sr 


, e = 


sr 


sr SG _ 


sr 


SG 


sr 


SG 


_ sr 


SL' iSN ~ 




iSP^ 


SP^' 


isT 


sT 



identities (46), (47) can be expressed in terms of the function F as 

sr sr 1 sr sr sr sr ^ sr 

SA SK 4 S^ SN Sr] SL St] ^ ' 

J^, sr sr sr sr , 
+ E(-<^'^> + <^'^» = 0' 
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Formulas (50), (51) take a slightly simpler form for the modified effective action f defined by 



T = T + ^<J^,J^> . 



Indeed, from (50), (51) we have 



ST 1 ST ST sr sr 

5A SK 4 S^ 5N Sr] SL ^ ' 

J^, Sf Sf ST ST , ^ 

^ Sr]i SP' Sr]i sp' 

^•^+^=°- 

For the operation in the r.h.s. of formula (52) we shall use the compact notation 

STi ST2 1 STi ST2 STi ST2 , , 

ri*r2=<— i,— 1>+-<-^,-|> + <-^,— 53 

SA SK 4 5$ SN Sr] SL 

, STi ST2 STi ST2 . 

.^^ Sr]i SP' Sr]i ^p' 

Now we can prove that the theory with generating function (8) is rcnormalizable. First recall 
that by simple dimensional counting the integrals contained in the perturbative loop expansion 

f = f ° + + f 2 + • . • (54) 

for the generating function T can be dimcnsionally regularized. Assume for a moment that 
this is done, i.e. the Wick rotation is performed and all the integrals are taken over R'^ as in 
the one-loop case discussed in the previous section. When the dimension d of the space R'^ 
goes to 4 each term in expansion (54) acquires a divergent term corresponding to a pole in 
variable £ = 4 — dase— >0 (minimal subtraction). The divergent terms T^^y can be determined 
recursively from the identity _ „ „ 

reg R divi 

where T^^g is computed taking into account all lower-order counterterms in the perturbation 
expansion, and is the renormalized effective action at the n-th order. 
At the zeroth order we obviously have 

f° = -i < F,F > +^ < > - J < > - < r],M^ > + 

3 



-I- ^ < dAVi, dAVi > + < K,sA> - < L,sv > N,s<^ > + (55) 

i=l 

3 

+ > + <P\ srii >) = 
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1 1 7TI 

= --<F,F> +— < 0^$,$ > -— <<^,F> + <{K-r]-(p),sA> + 
+ ^ < dAVi, > - < L,sr]> +-< N,s^> + 

i=l 

3 

+ ^(< P\ sr]i> + < srJi >), 

i=l 

where rj ■ (j)^ = ^"(/'^f,- 

From (52) we infer by induction over n that every divergent part F^jj, must satisfy the following 
equation (see [5], Sect. 12.4.3 for details): 

f°*rs„ + r3„*fo = o. (56) 

Equation (56) allows to find the most general form of the divergent terms. 
If we assume, by the choice of gauge, that the symmetry of FJJj^ under constant gauge trans- 
formations is^ unbroken then simple power counting and the ghost number conservation restrict 
the form of T^^^ to 

V^^^,=C^ + a<{K-r]-4>),8A>+{a-p)<{K-r]-cl,),[A,ri\>+ (57) 

1 ^ , 

+7 < L,ST]> +6- < TV, s$ > + + Y^isi < P\ srji > +ei < P\sr]i >), 

i=l 

where a, /3, 7, 5, Si, Si are numerical coefficients, and £" is a functional of degree four depending 

on A^<^,r|i,ri^ and having zero ghost numbers. 
Substituting expression (57) into (56) we get 

13 = 5 = £i=Si = -7, (58) 



5C SC^ 1 

(a -/3) < > + < >+-<—,.$> - (59) 

- < -T—,sr]i > - < -T^,sVi >= 0, 
6r]i 6r]i 

where 

11 ^ 
C = --<F,F>+—<nA^,^>-^<^,F>+J2<dAr]i,dAVi>- 

i=l 

The general solution to equation (59) is 

N / . 1 SC SC 5jC ^ , ^ 

= {P-a){< —,A> +-<_,$>_< > - < —fj, >) + (60) 

whore dnv is an arbitrary gauge invariant functional of degree four depending on A, $,7?8,r?j 
and having zero ghost numbers. 

We actually do not need to consider the most general expression for Cinv Some terms in this 
expression are a priori not required for the renormalization. From formula (7) it follows that for 
ni = the spatial components of the ghost tensor field $ do not interact with each other and 
with the other ghost fields. Taking into account the ghost numbers conservation we infer that 
the counterterms containing interactions of different components of $ and of anticommuting 
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ghosts r]i, rji for different indexes i should appear with multiples m or ^. The latter possibility is 

not realized in case of the minimal subtraction. Note also that the ghosts rji, rj^ symmetrically 
appear in formula (7). Therefore they must symmetrically appear in the expression for the 
counterterms. 

These rules restrict the expression for required for the renormalization to 

A™ = -y < ^'^ > +i < > < > + (61) 

3 3 ^2 

+ ^ 02 < dAVi, dAVi > + "4m^ < rJi, -rji > +05^ < i>, $ >, 

i=l i=l 

where 00,01,02, 03, 04,05 are nimerical coefficients, and four new mass parameters m^, rrin 

arc, of course, proportional to m. 

Substituting (60) and (61) into (57), using relations (58) and the homogeneity property of 
the functional with respect to coupling constants and fields we finally obtain the following 
expression for the most general form of the counterterms (compare with [5], Sect 12.4.3., formula 

12-165) 

r3.. = {(/3-«+y)<A^>+4(/3-« + y)<$,^>+(a3-^-)m— + 

3 

+ ^(/3 - a + ^)(< > + < ^ >) + -12- 2(/3 - a))^^ + 



=1 

+(05 - Ol - 2(/J - " ^05^ + (/? - 2a + ^) < i, A >}(?«)' 



where 



-11 ^ 



32 ' 4 

3 2 

- < ^, 7W?7 > + X TO^ < ^j, r?i > +^ <$,$> + < if, sA > - < L, sr/ > + (63) 

1 ^ -• 

+ -<N,s^>+ Y,{< P\ sr]i> + < P\ srJi >). 

i=l 

Formula (62) tells that the divergent part of the effective action at the n-th order has the 
same form as action (63) with renornialized coupling constants, masses and fields. Comparing 
formula (63) with (55) we see that a priori new counterterms have to be added to the original 
action (55) to cancel all the divergences. The only difference between expressions (63) and (55) 
is that the ghost fields acquire masses. 

The renornialized generating function Gr{J) is obtained by taking the inverse Legendre trans- 
form of tlie renormalized effective action r/j(A, $, rj, rj, rji^rj^^ K, L, N, P\ P' ) and by putting all 
the arguements, except for J, equal to zero. Recalling the definition (48), (49) of the Legendre 
transform we get 

GR{JJ,^,le,t,K,L,N,P\T)=irR{A,^,r],r],Vi,%,K,L,N,P\T) + 

1 - ^ -i 

+i{< J,A> --< I,^> - <^,r]> - <^,fj >-Y,i< t,m > + < >)), 
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Gr{J) = Gr{J, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0). 
This completes the perturbative study of the theory with generating function (8). 
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